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Abstract 

We report on improved isotope shift measurements of the isotopes 40,42,44,48 Ca in the 4s 2 Sl/2 —> 4p 2 P 3 /2 
(D2) transition using collinear laser spectroscopy. Accurately known isotope shifts in the 4s Sl/2 ~> 4 p P1/2 
(Dl) transition were used to calibrate the ion beam energy with an uncertainty of A U « d=0.25 V. The accuracy 
in the D2 transition was improved by a factor of 5 — 10. A King-plot analysis of the two transitions revealed 
that the held shift factor in the D2 line is about 1.8(13) % larger than in the Dl transition which is ascribed 
to relativistic contributions of the 4 p 1 / 2 wave function. 


1 Introduction 

Hyperhne structure and isotope shifts have been stud¬ 
ied in various transitions of calcium isotopes for 
decades both experimentally (see e.g. l|[T2|) and the¬ 
oretically (e.g. 113f[l8]). This is due to the importance 
of this element in various fields: With 40 Ca and 48 Ca 
the calcium isotopic chain is the only one that con¬ 
tains two naturally occurring doubly magic isotopes. 
Moreover, the nuclear charge radii along these stable 
isotopes show a very peculiar behaviour, having a max¬ 
imum at 44 Ca, a strong odd-even staggering and almost 
identical charge radii for 40 Ca and 48 Ca [l9,20] despite 
the large difference in nuclear mass. The isotopic chain 
has therefore been of great interest for nuclear struc¬ 
ture physics. Additionally, the long-lived isotope 41 Ca 
can be used for dating applications 21 as well as a 
tracer in biology and medicine to study the physiology 
of calcium in the human body 
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applications in ion traps and has thus also been studied 
with high accuracy 10 28 29 . Moreover, the isotope 
shift information in the 3d 2 Dj -± 4 p 2 Pj infrared 
triplet 30 has led to the discovery of an anomalous 


isotopic composition in mercury-manganese (HgMn) 
stars, in which the isotopic Ca ratio in the stellar at¬ 
mosphere is dominated by 48 Ca [31 32 . 


Precision data are available for many transitions, but 
the isotope shift in the 4 s 2 Si/ 2 ~± 4p 2 P 3 / 2 transition of 
calcium, often called D2 linq^J is only known with mod¬ 
erate accuracy 33 , even though this transition is pre- 


Hence, isotope 
shifts and hyperfine structure were also studied in sev¬ 
eral atomic transitions 25p 6 for isotope selective reso¬ 
nance ionization and ultra-trace analysis [27 . The cal¬ 
cium ion is a workhorse in the field of quantum-optical 


ferred for studies of nuclear moments in short-lived iso¬ 
topes. It allows for an extraction of the nuclear spectro¬ 
scopic quadrupole moment from the electric hyperfine 
structure, in contrast to the 4 s 2 Si/ 2 4p 2 Pi/ 2 (Dl) 

transition. This is of special importance because of the 
revived interest of nuclear structure physics in spectro¬ 
scopic information on stable and short-lived calcium 
isotopes. The interest concentrates on the neutron- 
rich region where new shell closures have been observed 


1 since it is isoelectronic to the corresponding Fraunhofer D2 
line in Na 
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recently for the isotopes 52,54 Ca at neutron numbers 
N = 32 [34] and TV = 34 35l. These shell closures 


line’ collinear laser spectroscopy measurements. 


are associated with three-nucleon forces and can be 
studied theoretically using chiral effective field the¬ 
ory p7p6 . In order to provide additional information 
about these shell closures, nuclear moments and iso¬ 
tope shifts of neutron-rich calcium isotopes have been 
studied at ISOLDE/CERN. The COLLAPS collabo¬ 
ration has measured hyperfine structure and isotope 
shifts beyond N = 28, namely the isotopes 49-52 Ca. 
Results for the nuclear moments have already been re¬ 
ported 37 , while the analysis of the isotope shifts is 
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2 Experimental setup 

The TRIGA-LASER experiment is part of the TRIGA 
SPEC setup at the TRIGA research reactor Mainz (40 



hampered by the large uncertainties in the D2 isotope 
shifts of the stable isotopes. These data are required 
because of the large measuring uncertainty of the exact 
high voltage potential of the RFQ cooler and buncher 
ISCOOL 38||39 . This device is used to collect the ions 
delivered from the ISOLDE on-line ion source and to 
send them in short bunches towards the collinear laser 
spectroscopy setup. The ion beam energy is usually the 
limiting systematic uncertainty for isotope shift mea¬ 
surements in collinear laser spectroscopy. But if ac¬ 
curate isotope shifts are available for several isotopes, 
the beam energy can be calibrated using this informa¬ 
tion. In order to provide the required reference data 
for the analysis of the short-lived isotopes, we have per¬ 
formed isotope shift measurements of 40 > 42 > 44 > 4 8Q a + a ^ 
the TRIGA-LASER setup in the Institute of Nuclear 
Chemistry in Mainz 40 . High-accuracy data of the 
isotope shifts in the 4 si/ 2 4py 2 transition, recently 
obtained using logical spectroscopy in an ion trap with 
an accuracy better than 100 kHz 


were used for a 
reliable calibration of the ion source potential. In this 
way we have extracted the isotope shifts in the D2 line 
with up to tenfold improved accuracy. In addition, the 
two transitions D1 and D2 can be compared via a King 
plot to examine the effect of the electron density of the 
pi/ 2 electron at the nucleus. Similar studies have pre¬ 
viously been performed for instance in Ba 41 , Mg 42 
and Ra 
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Figure 1: Schematic drawing of the main parts of the 
TRIGA-LASER beamline. 1. ” high voltage cage” with 
offline ion source; 2. ion steering elements; 3. charge 
exchange cell (not used); 4. mirror system for fluo¬ 
rescence detection; 5. Photomultiplier tubes; 6. laser 
beam 

The measurements on stable calcium isotopes pre¬ 
sented in this work were performed on ion beams of 
typically 0.8 nA with the natural isotopic composition 
produced in an offline ion source: A graphite tube with 
a central bore was filled with a few grains of metal¬ 
lic calcium. This crucible was heated by applying a 
current of about 50 A causing a voltage drop of about 
4 V along the crucible. The complete ion source is 
mounted on a high-voltage platform to which a poten¬ 
tial of about 10 kV is applied. The ions are accelerated 
towards ground potential. The ion source voltage is 
generated by a Heinzinger PNChp60000-10ump power 
supply with a relative stability specified to be better 
than 10 -5 . The potential of the high-voltage platform 
is continuously monitored and recorded using a Julie 
Research KV-10R high-voltage divider. The divider 
ratio is known with a relative precision of 5 x 10 -5 and 
a high-precision voltmeter of the type Agilent 34401A 
is used to measure the voltage across the precision re¬ 
sistor at the end of the divider chain. The divider ra¬ 
tio has been determined by comparison with the KA- 


The reactor will provide exotic nuclei produced by 
neutron-induced fission inside the reactor core for ’on- 


TRIN high-voltage divider 44 45 at the University of 
Munster. The starting potential of the ions inside the 
graphite furnace can differ from the platform potential 
by up to 4 V due to the voltage drop caused by direct 
heating of the graphite furnace. A 10° deflector is used 
to superimpose the uv-laser beam with the ion beam in 
anticollinear geometry. The beamline is equipped with 
two pairs of deflector electrodes in vertical and hor¬ 
izontal direction and a quadrupole doublet for beam 
steering and shaping. Superposition of the two beams 
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is ensured using several apertures with different sizes 
at two positions of the beamline. 

The optical detection region is supplied with two mir¬ 
ror systems for light collection and appropriate photo¬ 
multiplier tubes. They are mounted on an insulated 
platform to allow for Doppler-tuning by accelerating 
or decelerating the ions into the detection region. An 
offset potential of up to ±10 kV can be applied to the 
platform by a high-voltage device with similar spec¬ 
ifications as described above. This supply is again 
mounted on a platform and can be raised in poten¬ 
tial by a Kepco BOP 500M fast high-voltage amplifier. 
This allows fast Doppler scanning across ±500 V, corre¬ 
sponding to almost 30 GHz scanning range in the case 
of calcium and 10 kV beam energy. The voltages ap¬ 
plied to this platform during the measurements were 
also continuously monitored and recorded by a high- 
voltage divider and an Agilent 34401 A. 

A Matisse TS titaniurmsapphire (Ti:Sa) ring-laser 
pumped by a Verdi Nd:YV04 laser is used to pro¬ 
duce infrared light at 794 nm and 786 nm for exci¬ 
tation of the 4si/ 2 -A 4py 2 and the 4s y 2 —>• 4 p 3 / 2 tran¬ 
sition, respectively, after second-harmonic generation 
in a Wavetrain frequency doubling cavity. The Ti:Sa 
laser is actively stabilized to an external cavity to cor¬ 
rect for short-term fluctuations and to a FC1500 fre¬ 
quency comb to ensure long-term stability and an accu¬ 
rately known laser frequency during the measurements. 
The linewidth of the laser is approximately 250 kHz per 
4 hours. The uv light is coupled to a 180 m long single¬ 
mode fiber and transported to the beamline. The beam 
leaving the fiber is collimated, beam-halo is removed 
by a spatial mode filter to reduce scattered light back¬ 
ground and coupled into the beamline through a quartz 
window. The light leaves the vacuum system through 
a second window on the opposite side of the beamline. 
The light power was kept below 150 /iW in a beam 
of approximately 1.5 mm diameter and no appreciable 
sign of saturation broadening was observed at these in¬ 
tensities. With the optical detection setup composed 
of two mirror systems based on non-imaging optics an 
efficiency of one detected fluorescence photon per 200 
calcium ions has been demonstrated for the even iso¬ 
topes without hyperfine splitting. 


3 Experimental procedure 

We have performed measurements in the 4s y 2 —» 4py 2 
transition before and after the measurements in the 
4sy 2 -> 4 p 3 / 2 line. All measurements were performed at 
a constant starting potential without changing the ion 
optical settings and with the same laser system. For 
switching between the two transitions, the laser was 
changed in wavelength and the SHG cavity was opti¬ 
mized again. The fiber transport of the laser beam en¬ 
sured that the alignment of the laser through the beam- 
line did not change. The laser frequencies for both 
transitions were chosen such that comparable voltages 
had to be applied to the optical detection region for 
each isotope in order to avoid systematic changes of 
beam focussing. To increase the measurement accu¬ 
racy, spectra of the heavier calcium isotopes were al¬ 
ways alternated with reference measurements of 40 Ca + . 
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Figure 2: Typical spectra of 40 Ca + and 48 Ca + . The 
upper panel shows the 4sy 2 —>> 4p 3 / 2 transition (D2) 
and the lower one the 4sy 2 4py 2 transition (Dl). 


4 Analysis 

Line profiles of 40 Ca + and 48 Ca + in both transitions 
are depicted in Fig.[2j In the analysis of the spectra, 
for each data point the total acceleration voltage was 
converted into the corresponding Doppler-shifted fre¬ 
quency in the rest frame of the Ca + ion. All lines were 
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then fitted with Voigt profiles using a fixed Lorentzian 
width of Az/l = 23.4 MHz (natural linewidth). Gaus¬ 
sian linewidth Az/q, amplitude, background and peak 
center position z/q were used as free parameters. The 
linewidth parameter Az /q ^ 17 — 20 MHz was obtained 
from the fitting routine, with little changes between 
the isotopes due to the different masses. 

A slight asymmetry that is observed in the peak shape 
was similar for all isotopes. Tests with different line- 
shapes that better fitted to the profile, did not result 
in any significant change in the isotope shift. 


4.1 Beam energy calibration 

To increase the measurement accuracy, resonance scans 
of the isotope of interest x Ca + (X = 42,44,48) and 
the reference isotope 40 Ca + were alternated such that 
two scans of x Ca + were interleaved by four scans of 
40 Ca + to reduce the influence of long-term fluctua¬ 
tions. The whole measuring cycle (cycle I) was re¬ 
peated at a day about three months later (cycle II). In 
the meantime the beamline was modified at different 
points and everything was readjusted. As described 
above, the starting potential is known to U sourC e-iy, 
where the uncertainties are determined by the accuracy 
of the voltage measurement (AU/U ~ 10 -4 ) and the 
voltage drop along the crucible that will always lead 
to a higher voltage than that applied to the platform. 
To reduce the uncertainty, the resulting isotope shifts 
5 z^jg’ 40 ( U start) of all isotopes in the D1 line for various 
starting potentials U start between 9996 V and 10 kV 
were derived and for each voltage the mean square de¬ 
viation from the precise calibration data 


X (Ustart) 


E 

A= 42,44,48 



(E4tart) 

as4 A0 
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(1) 

with the purely statistical uncertainty A£zyg , was 
calculated. The value of y 2 obtained as a function of 
the starting potential is plotted in Fig. [ 3 ] for the case 
of the second cycle. The y 2 data are very well fitted 
by a parabola and the minimum of the curve is the 
most probable starting potential. The uncertainty is 
estimated to be 0.25 V, which corresponds to an in¬ 
crease of y 2 by 1 as indicated in Fig. [ 3 ] The results 
are 9997.85(25) V for the first and 9998.10(25) V for 
the second measuring cycle. The small change in the 


starting potential is within the uncertainty of the re¬ 
sults and might actually be due to a shifted ionization 
region for the ions within the graphite furnace caused 
by, e.g., different filling levels with calcium. 

The results of the calibration measurements are pre¬ 
sented in the lower panel of Fig.[3J Shown are the 
statistical uncertainties. The absolute values of the 
isotope shifts can be found in Table [lj 
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Figure 3: In the upper panel the y 2 -values according to 
Eq.[l]are plotted for different assumptions of the accel¬ 
eration voltage. The deviation between the obtained 
isotope shifts in the D1 transition and the calibration 
data 29] (black squares, set to zero) are plotted for 
the different isotopes in the lower panel. Error bars 
are purely statistical uncertainties. 


4.2 Results 

With the calibrated starting potential we can derive 
the isotope shifts in the 4 si/ 2 —)> 4p 3 / 2 transition. 
All measurements of 40 > 42 > 44 > 48 Ca + of both cycles 
have been combined to yield the isotope shift values 
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listed in Table Q] The first bracket for each value 
shows the statistical uncertainty which is the standard 
deviation of all measurements and the second bracket 
represents the systematic uncertainty. The latter 
depends mainly on the two voltage uncertainties: the 
one of the starting potential Af7 s tart and the one of 
the voltage floating the optical detection region Af7 0 d- 


The systematic uncertainty is derived as 46 


of the respective acceleration voltage of AU od /U od = 
5 x 1(T 5 . 

The determined isotope shifts are shown in Fig.[4]where 
the error bars represent the total uncertainty which 
was obtained as the geometrical average of the statis¬ 
tical and the systematic uncertainty. The agreement 
with literature as well as the significant improvement 
in accuracy is clearly visible. 


A5v, 


A,A' 

IS 


vlMA^^-AU + AM), (2) 


with 


au = 1( it l 

^ start 


5m 

m 


2 me 2 

Ab^start 

bis tart 


U 0 d Ai/od 


U R i 


u od 


and 


AM = 


AmA Am\' 


(3) 
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m m 

and the reference mass m, the mass difference 5m = 
77iA 7 — 777 a, the speed of light c and the laser frequency 
uy. The isotope masses m A A' have been derived from 
the AME2012 atomic-mass evaluation table 47 , tak¬ 


ing into account the mass of the missing electron while 
the missing binding energy of the electron needs not to 
be considered here. 
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Figure 4: Deviations of D1 and D2 isotope shifts mea¬ 


sured in this work and in Martensson-Pendrill et al. 
from the present values trig A- 
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Table 1: Isotope shift measurements in two transitions 
of singly ionized calcium with respect to 40 Ca + . The 
first bracket gives the statistical, the second one the 
systematic uncertainty. 


A 

< ^'4si /2 ->-4pi /2 

(MHz) 

< ^* y 4si /2 -!-4p3 / 2 

(MHz) 

Ref 

40 

0 

0 


42 

427.0 (11)(13) 

425(4)(4) 

425.706(94) 

426.4(15)(10) 
430 (18)(4) 
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29 

44 

848.8(9) (18) 

842(3)(8) 
849.534(74) 

850.1(10)(19) 

842(13)(8) 
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S3 

48 

1707.2 (29) (38) 

1696(4)(16) 

1705.389(60) 

1710.6(35)(39) 

1699(19)(16) 
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A A! 

The isotope shift 5 zq s ’ between two isotopes A and 
A' is caused by the change in mass and volume of the 
nuclei and can be split into the mass shift (MS) and 
the field shift (FS). 


The functional relation 
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fys — + 5v v 
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= F ■ 6 (r 2 ) A,A ’ + M ■ 
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rriA' — 777/A 

m A • m A ' 


(5) 


(6) 


with field shift factor F, mass shift constant M and 

change in mean square charge radius 5 (r 2 ) ’ is well 
known (see e.g. [33]). Both, the field and the mass shift 

depend on the investigated transition, but 5 (r 2 ) 
is identical in both transitions. Hence it can be elim¬ 
inated by comparing two transitions i and j. Eq. © 
can thus be transformed into a linear relation between 
the isotope shifts in two transitions [48 49 : 


A A ’ A A ’ Fi 
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( 7 ) 


The offset voltages were measured continuously us¬ 
ing a voltage divider providing a relative uncertainty 
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8v 


A,A' 

i/j 


m A • m A ' x A A ' 

- dv. . 

m A - m A > 


( 8 ) 


is the modified isotope shift. 

The modified isotope shifts in the D2 line determined 
here are plotted in Fig. [5] as a function of the mod¬ 
ified isotope shift in the D1 line as measured in the 
ion trap 29 . The data show an excellent linear be¬ 


haviour as expected from Eq. [7] From the linear re¬ 
gression we obtained the ratio of the field shift factors 
of Fd 2 /Fdi = 1.018(13) which deviates from unity by 
1.5 a. The field shift factor F is proportional to the 
difference of the electron probability density at the nu¬ 
cleus 


l ^| 2 = 




4s 


1 /2 




4 P 1 /2 3 /2 


( 9 ) 


and might therefore be slightly smaller for D1 since the 
pi/ 2 electron can obtain a small admixture of s y 2 to 
the wavefunction. This would result in a ratio slightly 
larger than 1 as observed. For a firm statement, a 
higher precision is needed. 


as well compatible with zero within 1.3cr. In order to 
improve the accuracy of the field shift factor ratio and 
the mass shifts, more precise measurements in a Paul 
trap will be carried out in the D2 line 
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5 Conclusion 

We have performed precision measurements of the iso¬ 
tope shifts in the D1 and the D2 transition in singly 
ionized calcium. The results agree with existing lit¬ 
erature values and improve the accuracy in the D2 
line considerably. They provide important reference 
data for isotope shift values of stable calcium ions 
and helpful calibration data for collinear laser spec¬ 
troscopy measurements on the neutron-rich isotopes 
beyond 48 Ca investigating the 4 s 2 Pi / 2 -A 4 p 2 P 3 / 2 tran¬ 
sition. 

The field shift in the D2 line is 1.8(13) % larger than in 
the D1 line. This observation might be used as a bench¬ 
mark for relativistic atomic structure calculations. 



<K,di (GHz u) 


Figure 5: King plot showing the linear relation of the 
modified isotope shifts S zqs in the 4sy 2 -A 4p 3 / 2 (D2) 
and the 4 Si/ 2 -A 4py 2 (Dl) transitions. The insets en¬ 
large the data points to show the agreement with the 
fitted line. 

From the y -axis intersection we can obtain the 
” weighted” difference between the specific mass shifts, 
which is 

F D O 

M D2 - M m • = -6.4(51) GHz u (10) 

t Dl 
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